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Abstract. — A  multiple-regression  model  is  presented  that  -reflects 
statistically  significant  correlations  between  defoliation  by  the  gypsy 
moth,  the  dependent  variable,  and  a  series  of  biotic  and  physical 
independent  variables.  Both  possible  uses  and  shortcomings  of  this 
model  are  discussed. 


^HE  GYPSY  MOTH,  Porthetria  dispar 
(L. ),  is  one  of  the  most  destructive  pests 
of  hardwood  forests  in  the  northeastern  United 
States.  Populations  of  this  insect  have  some- 
times achieved  densities  of  more  than  5,000 
egg  masses  per  acre  and  have  maintained 
these  high  densities  for  several  years  running 
(Campbell  1973a).  Such  populations  are 
capable  of  killing  a  high  proportion  of  the 
overstory  trees  through  repeated  and  heavy 
defohations  [Kegg  1971). 

On  the  other  hand,  the  proportion  of  the 
foliage  that  will  be  removed  by  gypsy  moth 
larvae  in  any  given  place  and  year  depends 
not  only  on  the  number  of  insects  that  are 
present,  but  also  on  other  variables.  A  par- 
ticular insect  density  does  not  determine  a 
particular  defoliation  level.  Rather,  an  array 
of  defoliation  levels  may  result  from  that 
density  (Campbell  1966).  And  before  we  can 
predict  defoliation  for  any  given  time  and 
place  with  reasonable  accuracy,  it  is  clear  that 
relationships  between  defoliation  and  other 
heretofore  unidentified  or  at  least  unquanti- 
fied  variables  must  be  accounted  for. 


OBJECTIVE 

Our  objective  is  to  describe  a  multiple- 
regression  model  that  reflects  statistically  sig- 
nificant correlations  between  defoliation  by 
the  gypsy  moth,  the  dependent  variable,  and 
a  series  of  biotic  and  physical  independent 
variables.  Although  the  data  used  to  develop 
this  model  were  collected  from  an  area  in 
northeastern  New  England  between  1911  and 
1931,  it  seems  reasonable  to  suppose  that 
most  of  the  underlying  mechanisms  through 
which  a  particular  defoliation  regime  is  de- 
termined today  are  not  much  different  from 
those  in  1911-31.  On  the  other  hand,  the  sys- 
tem may  have  changed  since  that  time  to  the 
point  where  the  model  presented  is  of  little 
or  no  practical  immediate  value.  In  any  case, 
this  model  will  be  tested  for  predictive  accu- 
racy in  the  near  future. 

PROCEDURES 

The  biological  data  used  here  were  accu- 
mulated by  personnel  of  the  old  Melrose 
Highlands  Gypsy  Moth  Laboratory  in  Massa- 
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chusetts.  Procedures  used  to  process  these 
data,  and  a  complete  list  of  the  independent 
variables  tested,  have  been  described  (Camp- 
bell 1973b). 

Processed  data  included  an  estimate  of 
each  of  the  following  variables  for  each  plot 
in  the  Melrose  system  during  each  year  the 
plot  was  examined. 

Significant  independent  variables: 


ZDL 
ZDH 

MLSTYR 
MNOW 

Trend 
DOM 

TFOL 
DCA„ 

PM 


:  Percentage  of  the  plots  examined 
in  a  zone  in  the  fall  of  year  (n) 
that  contained  more  than  500 
egg  masses  per  acre. 

:  Percentage  of  the  plots  examined 
in  a  zone  in  the  fall  of  year  (n) 
that  contained  more  than  5,000 
egg  masses  per  acre. 

:Egg  masses  per  acre  in  the  fall  of 
year  (n— 1).* 

:  Egg  masses  per  acre  in  the  fall  of 
year  (n).* 
MNOW 


MLSTYR 

:  Percentage  of  the  dominant  trees 
in  the  stand  that  were  favored  as 
food  by  the  gypsy  moth  (FCA). 

:An  index  of  total  foliage  per  acre. 

:  Percentage  defoliation  of  favored 
food  trees  (FCA)  during  year 
(n). 

: Precipitation  in  May  of  year  (n), 
in  inches. 


Dependent  Variable: 

DEF  (n+i)  =  Percentage  defoliation  of  the 

whole  stand  during  year  (n  +  1). 

First,  the  above  data  were  stratified  to  rep- 
resent only  stands  where  more  than  one-half 
of  the  overstory  was  composed  of  oaks,  be- 
cause defoliation  in  this  stratum  is  of  great 
concern  to  many  people.  Other  strata  were  set 
aside  for  examination  later. 

Second,  the  dependent  variable,  defoliation 


(  "One  egg  mass  was  always  added  to  observed  egg- 
mass  density  per  acre  because  no  egg  masses  were 
observed  in  some  of  the  Melrose  plots  during  some 
years.) 


in  year  (n  +  1)  (DEF(n+i)),  was  transformed 
to  ln(101-DEF(n+i)).  Although  defoliation  in 
most  of  the  plots  ranged  between  zero  and  25 
percent  during  most  of  the  years  they  were 
observed,  this  insect  is  likely  to  damage  trees 
only  when  defoliation  is  greater  than  50  per- 
cent; and  the  above  transformation  seemed  a 
reasonable  way  to  increase  the  accuracy  of  the 
model  for  forecasting  these  higher  defoliation 
levels.  By  the  same  token,  of  course,  the 
model  is  probably  less  accurate  than  another 
transformation — for  example  ln(DEF(n+i))  — 
would  have  been  for  forecasting  situations 
where  little  or  no  defoliation  might  be  ex- 
pected. 

A  few  relationships  were  known  to  exist 
between  defoliation  and  the  independent  vari- 
ables. Others  were  merely  suspected.  Non- 
linear relationships  were  presumed  to  be  com- 
mon. Similarly,  interactions  were  suspected 
between  combinations  of  the  independent 
variables  and  subsequent  defoliation.  So  a 
variety  of  multiple-regression  models  were 
tested,  having  the  general  form  Y  =  b|^  +  b^Xj 
+  h,^X2+  ■  •  •  +'^nX^;  and  due  provision  was 
made  for  testing  for  both  non-linearity  and 
interactions. 


RESULTS 

The  equation  presented  in  table  1  repre- 
sents our  best  prediction  model  for  defolia- 
tion, using  the  criteria  of  least  squares. 
Statistical  relationships  between  defoliation 
and  independent  variables  cannot  be  visual- 
ized just  by  inspection  of  the  separate  terms 
in  the  model  because  both  the  dependent 
variable,  DEF(n+i),  and  current  egg-mass 
density,  MNOW,  were  transformed  to  In  (101 
-DEF,„+i,)  and  In(MNOW),  before  the 
model  was  developed.  On  the  other  hand,  the 
contribution  of  each  term  in  the  model  to  ex- 
pected ln(101— DEF,n+i)  )  can  be  examined 
separately,  and  the  expected  value  of  defolia- 
tion itself  can  be  determined  by  simply  trans- 
forming the  expected  value  of  In  (101— 
DEF,„+i))  back  to  DEF,n+i).  This  procedure 
has  been  followed  in  developing  figures  1  to  5. 
Since  these  figures  were  all  derived  by  using 
the  mean  value  of  each  independent  variable 
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Table  1. — Analysis  of  variance  in  defoliafion  in  year  (n  +  1 )  1  n  (1 01  —  DBF 
(n  +  D)  as  a  funciion  of  the  specified  independent  variables  (N  =  533) 


Intercept  —  4.5315 
Variable 

Regression 
coefficient 

Standard  error 
of  coefficient 

t  ratio  on 
524  d.f. 

(DOM)  •  (InMNOW) 

0.00098977 

0.00030319 

3.3* 

(DMO)  •  (lnMNOW)2 

-.00019212 

.00003268 

-5.9* 

(ZDL)  •  (ZDH)-' 

-.00000219 

.00000056 

-3.9* 

(DCA(n)  )  •  (PM)  •  (InMNOW) 

.00009470 

.00002394 

4.0* 

(TREND)  •  (InMNOW) 

-.05556634 

.01332184 

-4.2* 

(TREND)  •  (InMNOW) 2 

.01473679 

.00329108 

4.5* 

(TREND)  •  (lnMNOW)3 

-.00097398 

.00019814 

-4.9* 

(TFOL) 

.00091184 

.00034688 

2.6* 

*  Significant  at  0.01  probability  level. 


R2  =  0.41. 


that  was  not  otherwise  specified,  these  mean 
values  are  shown  below: 


Variable 
Low  zone  density  (ZDL) 
High  zone  density  (ZDH) 
Prior  egg-mass  density 

(MLSTYR) 
Current  egg-mass  density 
(MNOW) 


Mean  value 
53  percent 
14  percent 
2,705  egg  masses 

per  acre 
2,488  egg  masses 
per  acre 


Dominance  of  class-A  trees       86  percent 
(DOM) 

Index  of  total  foliage  61 
(TFOL) 

Prior  defoliation  (DCA(n))      22  percent 
Precipitation  in  May  of  year    3.17  inches 
(n)  (PM) 

Zone  densities  (ZDL  and  ZDH). — Defolia- 
tion in  any  particular  place  within  a  zone 
could  be  expected  to  increase  in  year  (n  +  1) 
as  the  percentage  of  the  acres  examined  within 
that  zone  in  the  fall  of  year  (n)  that  con- 
tained high  egg-mass  densities  increased.  This 
relationship  is  examined  in  figure  1  for  zones 
where  at  least  40  percent  of  the  acres  exam- 
ined in  the  fall  of  year  (n)  contained  more 
than  500  egg  masses  per  acre,  ZDL,  and  where 
at  least  10  percent  of  these  same  acres  con- 
tained more  than  5,000  egg  masses  per  acre, 
ZDH. 

Presumably,  the  relationships  in  figure  1  re- 
flect primarily  on  an  abundant  supply  of 
newly  hatched  larvae  within  the  zone,  and  on 
the  tendency  of  these  larvae  to  disperse  to 
other  locations.  In  addition,  Leonard  (1967) 
has  suggested  that  newly  hatched  larvae  are 
more  likely  to  disperse  when  they  are  crowded 
than  when  they  are  not,  and  has  indicated  a 
mechanism  through  which  this  might  be  de- 
termined. These  relationships  may  also  reflect 
a  tendency  by  mobile  predators — such  as 


flocking  birds — to  concentrate  on  high-level 
host  populations  when  they  are  available,  thus 
allowing  lower-density  populations  to  survive 
in  larger  numbers  than  they  otherwise  might, 
and  consequently  to  subject  the  stand  to 
more  defoliation. 

Egg-mass  densities  (MLSTYR  and 
MNOW). — The  relationship  between  current 
egg-mass  density,  MNOW,  and  defoliation, 
DEF,n+i,,  was  straightforward  as  long  as  the 
difference  in  density  between  last  year's  den- 
sity, MLSTYR,  and  the  current  level, 
MNOW,  was  not  great  (fig.  2).  Although  de- 
foliation clearly  increased  as  current  density 
increased,  the  figure  simply  does  not  support 
a  long-held  concept  of  gypsy  moth  popula- 
tion management,  which  is  that  heavy  defolia- 
tion can  be  expected  at  or  above  500  egg 
masses  per  acre.  Indeed,  even  when  cvirrent 
density  was  about  5,000  egg  masses  per  acre, 
the  model  indicates  that  only  about  38  per- 
cent defoliation  could  be  expected  if  the 
population  had  been  stable  at  5,000  egg 
masses  per  acre  for  2  years  running. 

Thus,  unless  the  model  is  greatly  in  error 
for  this  density  range,  either  most  egg-mass 
sui-veys  for  proposed  control  work  greatly 
underestimate  actual  current  egg-mass  den- 
sity, or  the  usual  decision  to  exercise  popula- 
tion control  at  500  egg  masses  per  acre  should 
be  re-evaluated. 

On  the  other  hand,  the  model  indicates  that 
egg-mass  populations  as  low  as  500  current 
egg  masses  per  acre  might  result  in  heavy 
defoliation  of  the  stand  following  explosive 
increases   in    population   densities  between 
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Figure  1. — Relationship  between  defoliation 
during  year  (n  +  1),  and  the  percentages  of 
acres  in  the  zone  in  the  fall  of  the  year  (n)  con- 
taining more  than  (a)  500  egg  masses  per 
acre  (ZDL),  and  (b)  5,000  egg  masses  per  acre 
(ZDH).  All  other  independent  variables  have 
been  held  constant  at  their  mean  values. 
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Figure  2. — Relationship  between  defoliation 
during  year  (n  +  1)  and  egg-mass  densities 
in  the  fall  of  year  (n)  and  (n  -  1),  (MNOW 
and  MLSTYR).  Outbreak  conditions. 
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Figure  3. — Relationship  between  defoliation 
during  year  (n  +  1),  egg-moss  density  in  the 
fall  of  the  year  (n),  (MNOW),  and  the  per- 
centage of  the  dominant  trees  in  the  stand 
that  were  in  food  class  A  (FCA).  Outbreak 
conditions. 
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Figure  4. — Relationships  between  defoliation 
during  year  (n  +  1),  egg-mass  density  In  the 
fall  of  year  (n)  (MNOW),  and  an  index  of 
total  foliage  per  acre  (TFOL).  Outbreak 
conditions. 
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years  (n— 1)  and  (n).  Since  many  control 
operations  have  been  aimed  at  just  such  popu- 
lations, it  seems  reasonable  to  suppose  that 
current  population  management  philosophy 
for  this  insect  may  have  been  dominated  by 
this  latter  fact.  In  any  case,  though  we  have 
observed  the  tendency  of  rapidly  increasing 
populations  to  cause  much  heavier  defoliation 
than  their  relatively  stable,  but  equally  dense, 
counterparts,  we  cannot  propose  a  specific 
mechanism  to  account  for  it.  Perhaps  it  is 
determined  primarily  by  a  time  lag  for  the 
establishment  and/or  build-up  of  the  natural 
enemies  of  the  insect;  or  it  may  reflect  pri- 
marily on  the  physiological  state  of  the  insect 
itself.  Possibly  mechanisms  of  both  sorts  are 
involved. 

Dominance  of  class  A  trees  (DOM). — Re- 
lationships between  the  percentage  of  the 
dominant  trees  in  the  overstory  that  were 
favored  as  food  by  gypsy  moth  larvae,  DOM, 
current  egg-mass  density,  MNOW,  and  sub- 
sequent defoliation,  DEF(n+i),  are  somewhat 
puzzling  at  first  glance  (fig.  3),  because  the 
dominance  of  class  A  trees  was  more  closely 
related  to  defohation  when  current  density 
was  high  than  when  it  was  low.  This  differ- 
ence may  reflect  the  imposition  of  a  food  limi- 
tation on  the  insect  as  density  increased. 

Foliage  index  (TFOL). — Defoliation  tended 
to  decrease  as  the  total  amount  of  foliage 
available  to  the  larvae  hatching  from  a  given 
egg-mass  density  increased  (fig.  4).  This  is 
only  logical  because  the  proportion  of  the 
total  foliage  consumed  per  larva  was  certain 
to  decrease  as  the  total  amount  of  that  foliage 
increased.  However,  defoliation  only  de- 
creased by  about  3  to  5  percent  as  the  index 
of  total  foliage  doubled,  which  may  indicate 
either  that  dispersal  losses  among  the  newly 
hatched  larvae  tended  to  be  greater  within 
poorly  stocked  stands,  or  that  much  more 
foliage  was  consumed  per  larva  as  the  total 
amount  of  available  foliage  increased,  or  both. 

Prior  defoliation  (DCA(„))  and  precipitation 
(PM) — Defoliation  in  year  (n+1)  was  re- 
duced as  a  function  of  prior  defoliation  of  the 
class  A  trees  during  the  preceding  summer, 
DCA,n);  and  more  defoliation  could  be  ex- 


pected in  year  (n  +  1)  when  precipitation  in 
May  of  year  (n)  was  low  than  when  it  was 
high  (fig.  5).  Probably  this  latter  relation- 
ship reflects  a  greater  tendency,  in  part,  for 
a  disease  epizootic  to  develop  in  the  insect 
population  issuing  in  year  (n  +  1)  as  a  con- 
sequence of  the  higher  humidity  conditions 
present  during  the  larval  stage  in  year  (n). 

In  fact,  neither  mechanisms  that  are 
associated  with  heavy  prior  defoliation,  such 
as  reduced  fecundity,  nor  those  that  may  be 
associated  with  high  May  precipitation  levels 
in  year  (n),  such  as  increased  disease  inci- 
dence in  year  (n+1),  offer  a  completely 
plausible  rationale  for  the  results  shown  in 
figure  5.  Another  rationale  might  suggest  that 
mechanisms  similar  to  those  described  by 
several  authors  {Wellington  1965;  Morris 
1969)  may  alter  the  physiology  of  the  insect 
in  generation  (n  +  1)  through  directly  influ- 
encing the  maternal  parent  during  generation 
(n).  This  possibility  should  be  explored. 

DISCUSSION 

These  results,  in  our  opinion,  emphasize 
two  immediate  needs: 

First,  the  model  accounted  for  only  41  per- 
cent of  the  variation  in  ln(101— DEF  ), 
the  variable  chosen  to  reflect  defoliation  in 
OAK  stands.  This  may  imply  that  we  either 
do  not  have,  or  have  not  used,  a  truly  rele- 
vant index  of  a  fundamentally  important — 
but  currently  unknown — process  in  determin- 
ing a  defoliation  level.  Thus,  although  our 
primary  purpose  was  to  produce  a  prediction 
model  for  defoliation  by  this  insect,  our  re- 
sults emphasize  the  need  for  specific  studies 
to  identify  and  understand  the  processes 
through  which  particular  levels  of  defoliation 
by  this  insect  are  determined. 

Second,  agencies  charged  with  the  manage- 
ment of  this  population  system  should  review 
their  criteria  for  control.  Unless  the  popula- 
tion system  is  now  behaving  quite  differently 
from  what  is  indicated  by  this  model,  cur- 
rent criteria  for  control  may  sometimes  lead 
to  the  treatment  of  populations  that  would 
have  caused  Httle  defohation  if  they  had  been 
left  untreated. 
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